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ABSTRACT 

We report detections of two candidate distant submillimeter galaxies (SMGs), MM J154506.4—344318 
and MM J154132.7—350320, which are discovered in the AzTEC/ASTE 1.1 mm survey toward the 
Lupus-I star-forming region. The two objects have 1.1 mm flux densities of 43.9 and 27.1 mJy, and 
have Herschel/SFIKE counterparts as well. The Submillimeter Array counterpart to the former SMG 
is identified at 890 /xm and 1.3 mm. Photometric redshift estimates using all available data from 
the mid-infrared to the radio suggest that the redshifts of the two SMGs are Zphoto — 4-5 and 3, 
respectively. Near-infrared objects are found very close to the SMGs and they are consistent with 
low-z ellipticals, suggesting that the high apparent luminosities can be attributed to gravitational 
magnification. The cumulative number counts at 5'i.imm ^ 25 mJy, combined with other two 1.1-mm 
brightest sources, are 0.70to!34 beg“^, which is consistent with a model prediction that accounts for 
flux magnification due to strong gravitational leasing. Unexpectedly, a z > 3 SMG and a Galactic 
dense starless core (e.g., a first hydrostatic core) could be similar in the mid-infrared to millimeter 
spectral energy distributions and spatial structures at least at > 1". This indicates that it is necessary 
to distinguish the two possibilities by means of broad band photometry from the optical to centimeter 
and spectroscopy to determine the redshift, when a compact object is identified toward Galactic 
star-forming regions. 
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1. INTRODUCTION 

In recent years, remarkable progress has been made in 
finding the brightest (S’gso/^m ^ 70 mJy or S'l.imm ^ 
30 mJy) populat ion of submillime t er (submm) galax - 
ies (SMGs, e.g., iBlain et ^ 120021 : iCasev et al.l l20l4l 
via square-degree scale surveys using far-infrared (FIR) 
to millimeter (mm) single dish telescopes in space and 
on the ground. Thanks to their apparent high lu¬ 
minosity, usually with the aid of gravitational mag¬ 
nification, the brightest SMGs offer a unique oppor¬ 
tunity to_Jnyes^ate_ghysicaI__gro 2 ertiesofthe_JSM 
e.g.j iHarris et al. I2010t H vison et al. I I20I0I: iCox et al.1 


201j|_ Danielsonet,^ l2nilL I2ni3l IScott et all 120111 


Valtchanov et al. 

Mi 

Combes et alJl2012l:IDecarli et al. 

20121: llouo et al.ll2012l 

Ludu et al.1 120121: IBothwell et al. 


20I3[ lOmont et al.l 1201 3fl . resolved star-forming activity 


_ _ -__ ___ _ng_ 

(e.g.jSwinbank et al.ll2010HNegrello et al.ll2010ilFu et al.l 

20121 ~13ll and gas dynamics (e.g., Riechers et al.ll2011l : 


Rawle et al.ll2013tlMessias et ai.l[2014l) at the peak of the 

star-formation history of galaxies. 

Furthermore, the brightest SMGs located in the high- 
redshift tail of the SMG redshift distribution (median 
z ~ 2-3, e.g., IChaoman et al.l l2005t lYun et al] 120121 : 
iSwinbank et al.l [201411 at z = 3-6 provide a severe 
challenge to theories of galaxy formation and evolution 
(jBaugh et al.ll2005l: iGranato et al.l[20041200^ . The sur¬ 
face density of the brightest SMGs holds integrated infor¬ 
mation on when galaxies undergo intense starburst and 
how frequently strong gravitational leasing occurs, on 
which galaxy formation models depend. But, the source 
counts of brightest SMGs in the high-z tail are highly un¬ 
certain because of poor statistics, and it is indeed hard 
to find them; only a single SMG with S'l.imm ^ 30 mJy 
is expected within ~ 1 -10 deg^ from limited knowl¬ 
edge of current studies (IScott et al.l 120121 : iVieira et al.l 
IMO). Large-area surveys using LTersc/ieZ/SPIRE (250, 
350, and 500 fim) have recently rev ealed strongly lensed 
galaxies bright at 250-500 u m ('e.g.. lNegrello et 311120101 : 
iGonzalez-Nuevo et al.ll201^ . These surveys bring about 
a great change in our understanding of the statistical 
properties of the SMG population. The SPIRE bands, 
however, preferentially select SMGs at modest redshifts 
(z < 3), and the study of SMGs in the high-z tail (z = 3- 
6) is still far from being complete. 

The mm to long-submm wavel ength cameras such 
as S CUBA. SCUBA-2 ('85 0 am. IHolland et ahl IT^ 
[Ml . M AMBO-2 (1 .2 r nm, [Krevsa et al.lll998ll . AzTEG 
(1.1 mm. I Wilson et al.1 l2008al) . and Laboca (870 /rm, 
iSiringQ et al.1 I2009H are complementary to those FIR 
to short-submm surveys; For example, observations at 
longer wavelengths with the MAMBO a nd AzTEC cam¬ 
eras a nd the South Pole Telescope (^SPT. ICarlstrom et al.l 
1201 If) may exploit a uniform selection function in red- 
shift space, allowing the effici ent detection of the bright¬ 
est SMGs out to z > 3 ('e.g.. iLestrade et al.ll2009|. 120101 : 


llkarashi et al.ll20Tll : IVieira et al.ll2013HWeifi et al.[i20H . 

We have used the AzTEC camera to carry out large- 
area surveys toward Galactic star-forming regions, which 
cover ~ 3() deg^ of the sky in total with la sensitivities 
of 5-30 mJy beam~^. From mm/submm number counts 


(iNegrello et al.ll201Cll: Vieira et al.llMl IScott et al.ll2Ml 


iTakekoshi et al.ll20H . several detections of ultra-bright 
(■S'l.imm ^ 30 mJy) extragalactic sources at cosmo¬ 



Figure 1. The 6^ x 6' far-infrared to mm images of MM-J1545 
(left) and MM-J1541 (right). The contours show the AzTEC 
i.l-mm image, starting from 2a with a separation of 2cr. The 
background images show the Herschel/SPIKE 250, 350, 500-/i.m 
pseudo-color images. 


logical distances are expected within the survey area. 
This is complementary to earlier attempts to search for 
extremely-bright SMGs and to constrain the brightest 
end of extragalactic number count s by exploiting subm m 
maps from Galactic surveys fe.g.. IBarnard et ^120041) . 

In this paper, we report the AzTEC detections and 
multi-wavelength analyses of two 1.1-mm bright sources, 
MM J154506.4-344318 and MM J154132.7-350320 
(hereafter MM-J1545 and MM-J1541, respectively; see 
Figure [T]). These sources are found toward the Lupus-I 
star-forming region, a local (z = 0) molecular cloud, and 
indeed the close proximity of MM-J1545 to the molecular 
cloud misled us to classify it as a starless core when it 
was initially identified with the AzTEC 1.1-mm camera. 
Multi-wavelength data collected by subsequent follow¬ 
up observations, however, strongly suggest that they are 
galaxies at cosmological distances as presented in this 
paper, illustrating the importance of multi-wavelength 
imaging and spectroscopy of such objects. 

The paper is organized as follows. Section |2| describes 
observations and archival data obtained toward MM- 
J1545 and MM-J1541. In Section |3l we present the re¬ 
sults from the observations and photometric redshift esti¬ 
mates. SectionjUdiscusses the gravitational leasing effect 
on both objects, the brightest end of the 1.1 mm number 
counts, and the FIR-to-mm colors of the mm-selected 
sources. In Section [H we also explore a possibility that 
MM-J1545 would be a Galactic starless core. Finally, 
Section [5] summarizes our conclusions. 

Throughout this paper, we assume a concordance cos¬ 
mology with rim = 0.3, riA = 0.7, and Hq = 70 
km s“^ Mpc“^. 

2. OBSERVATIONS AND DATA 
2.1. AzTEC Observations and Arehival Data 

The 1.1-m m data were taken w ith the AzTEC bolome¬ 
ter camera (IWilson et al.1 l^08aD installed on the Ata- 
cama Submillimeter Exp eriment 10 m telescope (ASTE, 
lEzawa et al.l 12004 1200^ located at Atacama (altitude 
of 4860 m) in the Chilean Andes, during August 2007 
to December 2008. The weather conditions during the 
runs were excellent; The typical 220-GHz zenith opacities 
were in the range of 0.01-0.08. AzTEC/ASTE provides 
an angular resolution of 28" in full-width at half maxi¬ 
mum (FWHM). The complete description of the Lupus-I 
starless core survey will be reported elsewhere (Tsuk- 
agoshi et al., in preparation). 
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Table 1 

Multiwavelength counterparts to MM-J1545 and MM-J1541. 


Instrument 

Band 

Flux density 

MM-J1545 MM-J1541 

Unit 

Ref 

VLA 

20 cm 


1.3 ± 0.5 

mJy 

1 

VLA 

6 cm 

66 ±5 



2 

ATCA 

7 mm 

210 ±35 



2 

NMA 

2.7 mm 

< 5.9 


mJy 

2 

SMA 

1.3 mm 

20.8± 1.9‘> 


mJy 

2 

AzTEC 

1.1 mm 

43.9 ±5.6^^ 

27.1 ±5.0“" 

mJy 

2, 3 

SMA 

890 fim 

69.7± 12.l‘> 


mJy 

2 

SPIRE 

500 fim 

134.9 ± 11 . 9 s 

150.5 ± 12.8S 

mJy 

4 


350 fira 

109.4 ± 11 . 4 s 

169.5 ± 11.8S 

mJy 

4 


250 fim 

40.9 ± 12 . 7 s 

144.0 ± 18.0S 

mJy 

4 

MIPS 

24 fim 

< 0.3 

1.51±0.3R 

mJy 

5 

IRAC 

5.8 fim 

< 4.5 


MJy 

5 


3.6 fim 

< 0.7 


MJy 

5 

MOIRCS 

2.15 fim 

< 2.7‘i 


MJy 

2 

WFCAM 

1.64 fim 

< 7.8‘i 


MJy 

2 


1.26 fim 

< 10.2‘i 


MJy 

2 

2MASS 

Ks 


< 122® 

MJy 

5 


H 


< 95® 

MJy 

5 


J 


< 60® 

MJy 

5 

DSS 

I 


< 93^ 

MJy 

5 


R 


< 14^ 

MJy 

5 


B 


< 2.5^ 

MJy 

5 


References. — (1) The NRAO VLA Sky Survey (NVSS), (2) This work, (3) Kawabe et ah, in 
preparation; (4) Herschel Science Archive (HSA); (5) NASA/IPAC Infrared Science Archive (IRSA). 


^ The flux density not corrected for the flux bias due to confusion noises. The flux may be deboosted 
by only < 10% if the number counts are as flat as dN/dS oc S~^‘^ at > 10 mjy. 

^ Flux density estimated from the visibility fitting. _ 

^ The flux density retrieved from the catalog of the C2D survey (lEvans et al.ll200^ . 120091) Data Release 
4, in which the object is identified asSSTc2d J154132.7—350320 

^ The 2cr limiting flux with a 370 aperture. _ 

® The limit may be inaccurate due to blending of the near-IR object (see ^ 13. H and Table 
^ 3a limiting flux of the DSS/SERC survey. 

s The flux density measured with the Herschel Interactive Processing Environment IHIPE. [Qttll20in) 
command SourceExtractorSussextractor. Source confusion is not accounted for in the error. 


JVLA 6 cm 


ATCA 7mm SMA 1.3mm SMA 09Ojum MIPS 24/im IRAC 5.0^m IRAC MOIRCS Ks 


WCAM J 



* ; 










Relative offset (arcsec) 



Figure 2. (top) The 15^^ x 15^^ postage stamp images of MM-J1545, centered at the SMA 890-/im position. From left to right we 
show images at radio to near-infrared wavelengths. The contours represent the SMA 890-/im image, which are drawn at 4cr and 8cr 
(Icj = 3.0 mJy bearn”^). For the interferometer images, we also show the beam sizes with hatched ellipses. An object seen in near-infrared 
images, referred to as J1545B, is closely associated with the SMA source, but is offsets by 079. (bottom) The 15^^ x 15^^ postage stamp 
images of MM-J1541, centered at the MIPS 24-/im position. The images at mid-infrared to optical wavelengths are shown. The overlaid 
contours show the 24-^m image, which clearly shows the systematic offset (3'') from a near-infrared and optical source, J1541B. 
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Table 2 

Photometry of optical/near-infrared objects J1545B and J1541B. 


Instrument 

Band 

Flux density 

J1545B J1541B 

Unit 

IRAC 

5.8 /im 

59.0 ± 39.2^ 


MJy 


3.6 /im 

49.8±5.2‘>' 


MJy 

WISE 

4.6 /im 

< 100® 

253 ± 19® 

MJy 


3.2 /im 

< 90® 

260 ± 12® 

MJy 

MOIRCS 

2.15 /im 

25.0 ±2.1 


MJy 

WFCAM 

1.64 /im 

27.5 ± 10.1 


MJy 


1.26 /im 

17.2 ±5.9 


MJy 

2MASS 

Rs 


231 ± 28® 

MJy 


H 


237 ± 28® 

MJy 


J 


286 ± 34® 

MJy 

DSS 

I 

< 93 

261 ± 44"^ 

MJy 


R 

< 14 

275 ± 33‘i 

MJy 


B 

< 2.5 

182 ± 22"^ 

MJy 


^ From the catalog of the C2D survey llEvans et al.|[200a. [2009lj 
Data Release 4, in which the object is identified as SSTc2d 
J154506.4-344318. 

^ From the WISE All-Sky Source Catalog. 

From the 2MASS All-Sky Point Source Catalog, in which the 
object is identified as 2MASS J15413256—3503233. 

^ From the USNO-Bl.O Catalog UMonet et al.ll200^ j. The uncer¬ 
tainty includes the statistical and systematic photometric errors. 

® The limit is uncertain due to source confusion from an adjacent 
bright source. 

The redu ction procedure is d escribed by iScott et ^ 
(|2008li and iDownes et al] (|2012D . To extract point-like 
sources, the time-stream data were intensively cleaned 
using a principal component analysis (PCA) algorithm, 
and then mapped. The PCA cleaning works as a high- 
pass filter in the spatial frequency domain, and thus emis¬ 
sion extended significantly compared with the beam solid 
angle is fully filtered out. The FWHM of the point re¬ 
sponse function is 37" after an optimal filter is applied to 
surpress high spatial frequency noises. The pointing was 
checked every 1 hr using nearby radio quasars, resulting 
in the astrometric accuracy better than 3". Uranus and 
Neptune were used for flux density calibration, yielding 
an absolute accuracy better than 10%. The resulting 
root-mean-square (r.m.s.) noise over the mapped 4-deg^ 
region of the Lupus-I cloud is 5 mJy beam“^. 

We retrieve Level-2 images of i7ersc/ie//SPIRE from 
the Herschel Science Archive (HSA). The SPIRE data 
were obtained in 2011 January in parallel mode (Ob¬ 
servation ID: 1342213182) and are processed through a 
pipeline software with a standard processing generation 
(SPG) version of v8.2.1. The resulting 250-500 /rm im¬ 
ages cover a 2° X 2.3° region of Lupus-I with typical Icr 
noise levels o f 5-15 mJy beam~ ^. The details of the data 
are given bv lRygl et al.l (l2013fl . 

Eor near to mid-IR photometry, we refer to the pub¬ 
lic source catal ogs obtained wit h the Infrared Array 
Camera fIRAC. lFazio et ^1200411 and Multiband Imag- 
ing Photometer for Svitz er fMIPS. iRieke et all 120041 
onboard Spitzer, WISE (iWright et al.l " lOfl and the 
Two Micron All Sky Survey (2MASS, Skrutskie et al.l 
1200611 , which are available at the NASA Infrared Science 
Archive (IRSA). The Spitzer photometric data are pro¬ 
vided by the Spitzer Space Telescope “Erom Molecular 
Cores to Planet-fo r ming Disks” (C2D) Legacy Program 
(|Evans et al.ll200d . 1200911 . In addition, we also retrieve 


the basic calibrated data (BCD) from the Spitzer Her¬ 
itage Archive (SHA) in order to estimate the noise levels 
of the IRAC 3.6 and 5.8 pern and MIPS 24 /rm photome¬ 
try. The BCD of IRAC and MIPS are processed through 
masking, flat fielding, background matching, and mosaic- 
ing using a standard single frame pipeline on the Mopex 
software. The resulting r.m.s. noise levels at 3.6, 5.8 and 
24 /im are 0.2, 1.5 /rJy beam“^ and 0.1 mJy beam“^, 
respectively. The details on the Svitz e r obs ervations and 
catalogs are reported in lEvans et al.l (12009( 1. 

The NRAO Very Large Array Sky Survey (NVSS, 
(Condon et aLlllOQ^ i and the Digitized Sky Survey (DSS) 
data are also available at the positions of MM-J1545 and 
MM-J1541. 

2.2. Ancillary Data for MM-J1545 

For MM-J1545, w e used the Submillimeter Array 
(SMA. IHo et al.(l^004ll to measure its precise position and 
spatial extent, and then subsequently obtained multi¬ 
wavelength ancillary data that we describe below. 

2.2.1. The Submillimeter Array Observations 

MM-J1545 was observed with the SMA at 890 p,m in 
2010 January and at 1.3 mm in 2011 April. The 890- 
/rm observations were performed in the extended config¬ 
uration with eight antennas, which provided projected 
baselines ranging from 15 to 171 meters. The observing 
conditions were good (225-GHz zenith opacity of 0.05). 
The double side-band (DSB) receivers were tuned to a 
local oscillator (LO) frequency of 335.15 GHz, and the 
correlator provided 4-GHz band width in each sideband. 
The 1.3 mm observations were carried out in the compact 
configuration. The receivers were tuned so that ^^CO, 
^^CO and C^®0 J = 2-1 emission lines from the Lupus-I 
cloud (i.e., z = 0), as well as 1.3 mm continuum emis¬ 
sion of MM-J1545, can be imaged. The LO frequency 
was set to 224.86 GHz. The atmospheric transparency 
was again good. In both observing runs, two quasars 
J1626-298 (^ggo^m = 1.6 Jy) and J1454-377 (0.33 Jy) 
were used for complex gain calibration, while the pass- 
band response was calibrated using 3C273 and 3C279. 
The absolute flux density was scaled from the primary 
calibrator Titan. The accuracy of the flux calibration is 
estimated to be 15%. 

All of the data editing and calibration were performed 
using iDL-based standard routines in the MiR software 
package. The calibrated visibility data were imaged 
Fourier-transform ed) and deconvolved using MIRIAD 
Sault et'aT](1995ll tasks, invert and clean, respectively. 
In continuum imaging, the upper (USB) and lower side¬ 
band (LSB) data, eliminating channels where the local 
molecular lines are expected, were combined before imag¬ 
ing. The natural-weighted synthesized beam sizes at 
335 GHz and 225 GHz are l'.'49x 1'.'17 (the position angle, 
PA = -53?8) and 4"94 x 3"21 (PA = -5?7), respectively. 
The resulting r.m.s. noise levels at 335 and 225 GHz are 
3.0 and 0.76 mJy beam“^, respectively. 

2.2.2. Karl G. Jansky Very Large Array Observations 

The Karl G. Jansky Very Large Array (VLA) 6 cm 
data were taken in C-configuration in 2010 Novem¬ 
ber (project ID: lOC-226). The correlator was config¬ 
ured to provide 16 x 128 MHz subbands covering from 
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4.2 to 6.1 GHz. The data were calibrated using the 
VLA Calibration Pipeli neF^ which is based on the CASA 
data reduction package (jMcMullin et al.l[2007li . Imaging 
was also carried out using CASA employing the multi¬ 
frequency synthesis algorithm (spectral Taylor expan¬ 
sion) with nterms = 2 during the deconvolution, to take 
into account the spectral ind ex of the sources within the 
field (|Rau &: Cornwellll2011ll . along with Briggs weight¬ 
ing (robust = 0.5). The resulting synthesized beam is 
lOl'Sh X 3"22 at position angle 3?7, and the r.m.s. noise 
is 5.8 /rJy beam“^. 

2.2.3. The Australian Telescope Compact Array 
Observations 

The Australian Telescope Compact Array (ATCA) was 
used to take 7 mm continuum data in the H214 array con¬ 
figuration in 2013 October and in the H168 configuration 
in 2014 April (project ID: C2910). We used five tunings 
to fully cover the 33.4-50.4 GHz band to search for a red- 
shifted ^^CO (2-1) line (Taniguchi et ah, in preparation). 

The center frequencies were set to 35.25, 38.75, 42.20, 
45.40, 48.60 GHz, and two 2 GHz spectral windows of 
the CABB correlator were configured adjacently, result¬ 
ing in an instantaneous frequency coverage of 3.9 GHz 
for each tuning. The complex gain was monitored using 
a radio source B1541—375 (S'rmm = 1-1 Jy, 10?7 away 
from MM-J1545). Bandpass and delay calibrations were 
performed using 3C 279 once per night before the obser¬ 
vations started. Mars was used for absolute flux density 
calibration. The absolute flux density uncertainty is es¬ 
timated to be < 15%. The data were calibrated using 
MIRIAD and imaged using CASA with Briggs weighting 
(robust parameter of 0.5). We did not use the five longest 
baselines including the antenna at the W392 station be¬ 
cause of poor phase stability. The resulting synthesized 
beam size and r.m.s. noise level were 6'.'02 x 4'.'43 (PA = 
82?4) and 35 /xJy beam“^, respectively. 

2.2.4. Nobeyama Observations 

We also use the Nobeyama Millimeter Array (NMA) 
at Nobeyama Radio Observatory (NRO) to constrain the 
3 mm photometry. The observations were done during 
2010 April and May in the D configuration, where only 
five antennas were operational. The DSB receivers were 
tuned at 98.20 GHz (LSB, A = 3.05 r nm) and 110.20 
GHz (USB), and the UWBG correlator (jOkumura et al.l 
l2000fl with a 1-GHz bandwidth was used. We did not 
use the USB data because of poor qual ity. The visibil¬ 
ity da ta were calibrated with uvproc-ii (iTsutsumi et al.l 
I1997I1 . and then imaged using the AiPS (iGreisenI l2003ll 
task, imagr. The resulting beam size and r.m.s. noise 
level are 1672 x 974 (PA = —16?3) and 2.0 mJy beam“^, 
respectively. 

Since MM-J1545 is toward the edge region of a Galac¬ 
tic molecular cloud, it is necessary to carefully investi¬ 
gate whether the source is indeed extragalactic. We used 
the NRO 45-m telescope to observe Galactic ^^GO (1- 
0), ^^GO (1-0) and C^®0 (1-0) emission lines. The 
45 m observations were performed during 2010 January 
to April. The ^^GO and observations were carried 
out with the on-the-fly (OTF) mode of the multi-beam 

See https://science.nrao.edu/facilities/vla/data-processin| 


BEARS receiver (jSunada et al.l [2000ll and with the po- 
sition switching mode of the TlOO single-beam receiver 
(|Nakaiima et al.l[2008ll . respectively. The ^^GO map cov¬ 
ered a 17' X 17' region including MM-J1545 and a part of 
the Lupus-I cloud located north-east of MM-J1545. The 
^^GO OTF observations covered a 4' x 4' region cen¬ 
tered on MM-J1545. In both of the observations, the 
AG45 digital spectrometer was used. For the G^®0 ob¬ 
servations, we used the TlOO single-beam receiver and 
the acousto-optical spectrometers with a high-dispersion 
mode (AOS-H) in the standard position-switching mode, 
providing a spectral resolution of 0.10 km s“^ at 110 
GHz. Intensity calibration was done using the single¬ 
temperature chopper-wheel method, and the accuracy of 
intensity calibration is estimated to be 20 %. 

2.2.5. Near-infrared Observations 

We carried out near-infrared (NIR) imaging observa¬ 
tions of MM-J1545 with the Suba ru telescope equipped 
with the MOIRC S instrument (jlchikawa et alJ I2006I 
iSuzuki et al.l 1200811 on 2010 April 23. The observations 
were made with the ATg-band filter at A = 2.15 /tm, with 
a pixel scale of 0.12" pix“^. The total integration time 
was 2.9 ksec. The position and magnitude were cali¬ 
brated with several 2MASS point sources with Kg ^ 15 
within ~ 1'. The systematic uncertainties of the astrom¬ 
etry and magnitude are estimated to be O'.'l and 0.07 
mag, respectively. 

JH-band imaging observations we re performed with 
the Wide Field Gamera (WFCAM, iGasali et al.l l2007h 
attached to the United Kingdom Infrared Telescope 
(UKIRT) on Mauna Kea, during 2010 March 15- 
19. These observations were complemented with the 
MOIRGS Ks and Spitzer data. The sky was photomet¬ 
ric throughout these nights and seeing sizes were mostly 
077-171. A five-point dithering and a four-point micro¬ 
stepping were used in all observations to compensate for 
bad pixels and to recover full point spread function sam¬ 
pling. This results in proper sampling of the seeing size 
with the 074 WFCAM pixels. Each exposure time was 
10 sec, yielding total integration times of 2.8 ksec and 
2.4 ksec at J and H, respectively. All of the data were 
reduced in a standard manner. The astrometric uncer¬ 
tainties are less than O'.'l, and the photometric errors at 
J and H are 0.18 and 0.32 mag, respectively. 

3. RESULTS 
3.1. Individual Sources 

MM-J1545 and MM-J1541 are detected at 1.1 mm 
at signal-to-noise ratios (SNR) of 7.8 and 5.4, respec¬ 
tively. The flux densities of MM-J1545 and MM-J1541, 
which are measured for the PCA-cleaned map, are 43.9± 
5.6 mJy and 27.1 ± 5.0 mJy. We do not correct for possi¬ 
ble flux boosting due to underlying fainter sources since 
it is difficult to estimate the Bayesian prior from uncer¬ 
tain number counts at S'l.imm > 10 mJy. They are, how¬ 
ever, unlikely flux-boosted because such bright sou rces 
are extremely rare and, as we shall discuss in § 14.21 the 
shape of the number counts are appear to be much flat¬ 
ter than those at S'l.imm < 10 mJy, where flux boosting 
is significant. These two are amongst the four bright¬ 
est ever discovered in the AzTEC/ASTE campaign, one 
of which is a Si.imm = 3’ 7 mJy SMG at z = 3.39 re- 
ing/^ 0 Pf^gn^,y llkarashi et al.l (|20Hil (known as Orochi or 
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Projected baseline v'(u^+v^) (kX) 

Figure 3. The flux densities versus projected baseline length 
normalized by an observing wavelength (i.e., spatial frequency, 
y/u^ + for the 890-/^m SMA image of MM-J1545. The filled 
circles show the visibility amplitudes of the SMA counterpart. The 
amplitude is well fitted by a Gaussian with a point source (dash- 
dotted curve), suggesting an unresolved point- or cusp-like struc¬ 
ture embedded in a 2" extended component. Just for reference, 
we plot circular-symmetric Gaussian profiles with a total 890-/rm 
fiux of 69.7 mjy and FWHM of 075, l70, and 270 (dotted curves). 
We sho w the Fourier components o f a 0.04 Mq Bonnor—Ebert (BE) 
sphere IIEbertiri955l: [Honnorlll956l) as a realization of a prestellar 
core. The gas temperature and central H 2 density is assumed to be 
7.1 K and 5 x 10^ cm“^, respectively. We also show the predicted 
890-/im Fourier distribution of a first hydrostatic core with inclina¬ 
tion angle of z = 90° collapsed from a 0.3 Mq BE s phere, which is 
computed by radiation hydro dynamic simulations JTomida et al.l 
I201QI: fS^aigo fc TQmisakall2Qlll'l . while we conclude that MM-J1545 
is not likely a first hydrostatic core, bu t a 2 : 4-5 starburst galaxy, 

from multi-wavelength analysis. See ^ l4.4l for details. 

HXMMO 2 , iWardlow et alJ[2013D and the other is a 43- 
mJy source toward the peripheral field of the Small Mag - 
ellanic Cloud reported elsewhere (jTakekoshi et alJl2ni3[ l. 
Both of the sources toward the Lupus-I cloud are also de¬ 
tected in the SPIRE 250, 350, and 500-^ni bands, most 
of which have flux densities of >100-mJy. In this re¬ 
spect, both are very similar in FIR-to-mm brightness to 
the extre mely-luminous SMGs detect ed in the Herschel- 
ATL AS (IGonzalez-Nuevo et al.l[^012f) and HerMES sur¬ 
veys (|Wardlow et al.ll2013ll . The coordinates and results 
of photometry are listed in Table [TJ The postage stamp 
images are given in Figure [5] 

An important point that we must note is that the 
sources are found toward a Galactic {z = 0) molecular 
cloud. Multi-wavelength data require an extremely cold 
(< 10 K) and compact (<C 10") nature of the sources, 
which is too rare among Galactic star-forming objects 
and thus unlikely associated with the molecular clouds, 
but are located at cosmological distances. The spatial 
profiles are all consistent with a point-like source with the 
AzTEG (~ 37") and SPIRE (20-30") beams, unlike star¬ 
less cores found in Galactic molecular clouds, which are 
typically found to be ~ 0.1 pc correspondi ng to ^ 100" 
at th e distance to the Lupus-I cloud le.g.. lOnishi et ahl 
l2002f) . The SPIRE photometry places a constraint on 
the peak positions of dust spectral energy distributions 


(SEDs) at Aobs ^ 300-400 jxm, which indicates the effec¬ 
tive dust temperatures of Tdust/(l + -z) ^ 10 K. Given the 
dust temperatures found in starburst galaxies (> 40 K), 
this suggests that MM-J1545 and MM-J1541 are situated 
at z > 3. We will perform photometric redshift estimates 
in § 13.21 Below is a summary of the individual sources. 

3.1.1. MM-JI 545 

The SMA interferometric observations at 890 /rm and 

1.3 mm confirms the exact position at the J2000 equato¬ 
rial coordinate of (a, S) = (15^45“6!347, -34°43'18"18) 
with the uncertainties of ~ 0'.'09 and ~ O'.'IO at 890 /rm 
and 1.3 mm, respectively (see below for estimation of po¬ 
sitional uncertainties). The flux densities at 890 ^m and 

1.3 mm are 69.7 ± 12.1 mJy and 20.8 ± 1.9 mJy, respec¬ 
tively. The 890 /rm image is resolved with its 1" beam, 
and the beam-deconvolved source size fitted with a sin¬ 
gle 2-dimensional Gaussian using a miriad task uvfit 
is 1'.'2 X Of'Gd (PA = 37°). Figure [3] shows the visibil¬ 
ity amplitudes versus projected baseline length (i.e., the 
Fourier transform of a radial profile as a function of spa¬ 
tial frequency) for MM-J1545. The Fourier components 
are well expressed by a single Gaussian (271 ± Of'G in 
FWHM) with a constant offset (25 ± 6 mJy), implying 
a cusp-like compact structure embedded in a extended 
(r; 2") component. 

At the SMA position, VLA 6 cm and ATGA 7 mm 
emission is also detected at 66 ± 5 /rJy and 210 ± 35 /iJy, 
respectively, but it is not resolved with the VLA and 
ATGA beams. We do not detect 2.7 mm emission in 
the NMA image, which provides a ia upper limit of 5.9 
mJy. The counterpart was not detected in the MIPS 
24 ^m down to the 3cr limiting flux density of 0.3 mJy. 

In Figure HI we fit the SED from the infrared to the 
radio with a single-component modified blackbody (or 
graybody), KdB,y{v,T), where is the Planck func¬ 
tion and Kd = kqI'^ is the dust absorption coefficient 
which follows a power-law function of frequency i'. The 
FIR-to-mm part of the SED is well described with a sin¬ 
gle modified blackbody with an effective temperature of 
7dust/(l + z) = 7.1 ± 0.3 K and the emissivity index of 
P = 1.4±0.1. However, the 6-cm flux clearly exceeds the 
gray-body function. The spectral index over the 6 cm 
band is a = 0 ± 1 (the error is the Icr confidence inter¬ 
val), where Si, oc This is rather flat compared with 
the Rayleigh-Jeans slope of a = 2 -|- /3 = 3.4 ± 0.1, sug¬ 
gesting that the 6-cm flux arises from synchrotron and/or 
free-free emission. 

In the MOIRGS, WFGAM and IRAG images, a faint 
(ATs = 18.55 in the Vega magnitude system, S's.e^m = 
49.8 ±5.2 ^Jy) source is detected at 0'.'9 east of the SMA 
890-^m peak. Hereafter we refer to this NIR object as 
J1545B. Table [2] lists the optical to NIR photometry of 
J1545B. Figure [5] shows the centroids and uncertainties 
of multi-wavelength counterparts detected in the NIR 
(2.15 ^m, 3.6 ^m) and the submm to centimeter (890 ^m, 

1.3 mm, 7 mm, 6 cm). The positional uncertainty of 
each counterpart is estimated by adding statistical and 
systematic errors in quadrature. The statistical error 
is obtained from A6»stat ~ 0.6 6»beam/SNR, where 6»beam 
and SNR are the beam FWHM and signal-to-noise ratio. 
Given the small positional uncertainty of the SMA and 
MOIRGS (< 0.1"), this offset is significant and it is un¬ 
likely that J1545B is a counterpart to MM-J1545. The 
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Wavelength (/tim) Wavelength (/um) 

Figure 4. Radio-to-FIR photometric redshifts of (a) MM-J1545 and (b) MM-J1541. {top) The upper panels show the reduced-x^ as a 
function of redshift. The lower panels show the far-infrared (FIR) luminosity that minimize x^ at given redshift. To avoid contamination 
from a (possible) foreground lensing object, we only use photometric data at Aobs ^ 24 /xm. (bottom) The best-fit SEDs of the targets. 
The black circles show flux densities of the targets which used in the photo- 2 : estimates, while the open squares indicate those of nearby 
lens candidates, J1545B and J1541B, which are not used in photo- 2 : estimates. The dotted curves are the best-fit function of a modified 
black-body (gray-body). The measured flux densities of MM-J1545 are well fitted by a gray-body with Tdust/(1 z) = 7.1 K and /3 = 1.4, 
while the SEE) of MM-jl541 is described by a Tdust/(1 -\-z) = 8.5 K gray-body (/3 is fixed to 1.8 because of the small number of photometric 
data points). 


FWHM of the MOIRCS image of J1545B is the 0"97, 
whereas that of a nearby star (point spread function, 
PSF) is (y/69, suggesting that J1545B is intrinsically ex¬ 
tended (a PSF-deconvolved size of ~ Cf/GS) and likely a 
low-z galaxy along the line of sight toward MM-J1545. 
This close association of a foreground galaxy may grav¬ 
itationally magnify the background object, which can 
naturally explain the extreme flux density of MM-J1545. 
The possible gravitational lensing will be discussed in 
§[4T] 

Furthermore, we do not detect the J = 2-1 and 1- 
0 transitions of C^®0 toward MM-J1545, as shown in 
Figure [6l From this we put a meaningful constraint on 
a molecular mass of a possible Galactic dense gas core, 
suggesting that MM-J1545 is not of Galactic origin. The 
J = 2-1 and 1-0 transitions of G^®0 trace molecular gas 
with n(H 2 ) > 10^'^ and > 10^’^ cm“^, respectively, and 
universally seen associated with Galactic starless cores. 
The 3(7 upper limits on the main-beam temperature with 
a velocity resolution of 0.5 km s“^ are T^b < 0.3 K (G^®0 
J = 2-1) and Tmb < 0.07 K (C^®0 J = 1-0), yielding an 
upper limit to the dense gas mass under the local ther¬ 
modynamical equilibrium (LTE) of < O.OOSMq from the 


G^®0 (1-0) constraint. This LTE mass is much smaller 
than those found in Galactic starless cores. No com¬ 
pact nor emission is significantly detected at 

the SMA continuum position although ambient molec¬ 
ular gases are contaminated across the field of view of 
the 45-m and SMA maps. The absence of high-velocity 
components in spectra, which trace molecular out¬ 

flows from an accreting protostellar system, rules out any 
protostellar phases. 

The 6 cm emission is also critical to judge if the object 
is extragalactic; starless cores have neither synchrotron 
nor free-free emission unlike galaxies. The clear excess 
to the gray-body defined at 250 /rm to 7 mm and the 
rather flat spectral index in the 6 cm band excludes the 
possibility that the 6 cm signal is dominated by dust 
emission from a Galactic starless core. Thus, it is natural 
to suppose the object to be extragalacti c. W e will further 
discuss on the Galactic possibility in § 14.41 

3.1.2. MM-JI 54 I 

MM-J1541 is also unlikely a Galactic source, because it 
is well isolated from the main clouds {Ay > 2) of Lupus-I 
and meets the “off-cloud” criterion defined bv iRvgl et al.l 
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(|2013f) . It is significantly detected at 24 //m with Spitzer 
(SNR ft! 6 ) at (a, 5) = (15^41“32!706, -35°03'19"03), 
all consistent with a z ft 3 SMG, although no IRAC data 
are available. A small enhancement (2.5cr, 1.3 mJy) at 
20 cm is seen in the NVSS 20-cm image. As shown in Fig¬ 
ure m the FIR-to-mm part of the SED is well described 
by a single component gray-body with an effective dust 
temperature of Tdust/(1 + z) = 8.5 K, where we assume 
a dust emissivity index of /3 = 1.8 because of the limited 
number of photometric data points. 

A bright NIR/optical source is detected in 2MASS 
JHKs and DSS BRI images at {a, 5) = (15^41“32!56, 

—35°03'23'.'3), which is 3" south-west of the 24 /im cen¬ 
troid. We hereafter refer to this 2MASS/DSS object as 
J1541B. We also find a 3.2 and 4.6 ^m object in WISE 
data at this position. Unfortunately, we have no higher 
resolution images at 24 /xm or longer wavelengths, at 
which the emission likely comes from MM-J1541. But 
given that the 24 /xm source is detected at high SNR 
and is likely a counterpart to MM-J1541, we can place 
a constraint on its position; the statistical uncertainty 
is estimated to be 0'.'6 while the systematic error is 1'.'4, 
yielding a total positional uncertainty of 24 /xm image of 
approximately 1'.'5. The positions of J1541B measured 
in six 2MASS/DSS bands coincide with each other, and 
the astrometric accuracy is estimated to be better than 
(y.'3 for J1541EE3- Therefore, the offset between the 24 
/xm peak and the 2MASS/DSS position is significant at 
a ~ 2(7 level. Furthermore the shape of the SED is con¬ 
sistent with that of a low-z passive elliptical as discussed 
later in §[111 suggesting that J1541B may be a galaxy 
which lenses the background SMG MM-J1541 at z ^ 3. 

We list flux densities of J1541B in Tableland will fur¬ 
ther discuss the possible lensing in § 14.11 

3.2. Photometric Redshift Estimates 

We fit the SED models to the photometric data 
points at Aobs > 24 /xm to constrain their photomet¬ 
ric redshifts and FIR luminosities. We use SED tem¬ 
plates of well-chara cterized sta r burst galaxies; Arp 220, 

M 82, (GRASIL, iSilva et HI 119981 1. a composite of 
76 radio-identified SMG s with spectroscopic redshifts 
(|Michalowski et a h! I2010 D, and SMM J2135—0201 (the 
cosmic eyelash. ISwinbank et al.ll2010ll to search for mini¬ 
mum by simply varying the redshifts and luminosities 
of the SED templates. We consider 20% of an absolute 
flux density uncertainty for all photometric points in ad¬ 
dition to the statistical error. We also take i nto account 
the 1(7 confusion noise (|Nguven et al.ll20T^ for photo¬ 
metric errors of the SPIRE bands. 

Figure [4] shows the results of photometric redshift es¬ 
timates and the best-fit SEDs. The inferred redshift 
of MM-J1545 is z ft 4-5, although the derived photo¬ 
metric redshifts depends on the templates. Overall, the 
measured SED is in good agreement with the Arp 220 
and SMM J2135 templates although the fit with M 82 
is poorer than the others. The best-fit redshifts and 
68 % confidence intervals are z = 4.67^0 74 (Arp 220), 
5.66lJ;8^ (M 82), 4.06t(!;?? (average-SMG) and 4.201)5 

Explanatory Supplement to the 2MASS All Sky 
Data Release and Extended Mission Products, § 11. 2, 

http: //www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html 



A R.A. (arcsec) 

Figure 5. The Icr positional uncertainties of MM-J1545 mea¬ 
sured at near-IR (MOIRCS As, purple filled circle; IRAC 3.6/rm, 
light blue dashed circle) and submm to radio wavelengths (SMA, 
AzTEC, ATCA and VLA, elipses in red, black green, blue and or¬ 
ange). The uncertainties are estimated by adding statistical and 
systematic positional errors in quadrature. The background image 
and contours show the MOIRCS Ks and SMA 890 /i.m, respectively. 
The contours start from 2cr with an interval of 2cr and the negative 
flux densities are represented as dashed contours. The separation 
between the NIR and the SMA/JVLA sources is significant at > 2cr 
confidence level. 

(SMM J2135). The inferred FIR luminosities for the 
best-fit redshifts are then log(LFiR/T 0 ) = 14.3lo 2 
(Arp 220), 14.3to;2 (M 82), Id.l+g;^ (average-SMG), 

and 13.9lo 2 (SMM J2135). The dust temperature is 
Tdust — 36-40 K if z ft 4.1-4.7, which is similar to those 
found in SMGs le.g.. iKovacs et al]l2006ll . 

The redshift of MM-J1541 is estimated to be z ft 3, al¬ 
though the available photometric data are very limited. 
Among the templates, Arp 220 better reproduces the ac¬ 
tual SED than the others, and the 24 /xm detection is well 
accounted for by the 7.7-/xm feature of polycyclic aro¬ 
matic hydrocarbons (PAHs). The photometric redshifts 
obtained the templates are z = 2.96lg;4o (Arp 220), 
3.58t°;^^ (M 82), 2.65t°;^? (average-SMG) and 2.53t°;^® 
(SMM J2135). The inferred FIR luminosities for the 
best-fit redshifts are then log(LFiR/T 0 ) = 14.Ol));^ 
(Arp 220), 14.llg;5 (M 82), 13.9tR2 (average-SMG), and 
13.6lg'2 (SMM J2135). If we consider z ft 3.0, then the 
dust temperature is approximately 34 K, again consistent 
with those found in SMGs. 

4. DISCUSSIONS 

4.1. Possibility of Galaxy-Galaxy Lensing 

The extreme luminosities of MM-J1545 and MM-J1541 
are likely attributed to galaxy-galaxy lensing; In many 
cases, an elliptical galaxy seen in the optical to NIR is 
associated with brightest (Tfir > 10^‘^Lq) SMGs. Inter¬ 
ferometric imaging of these bright sources often reveals 
multiply split images or extended structures that are well 
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(a) ^^CO(l —0) with 1.1mm 
(K km s“^) 
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(b) ^®C0(l-0) and C^®0(1-0) spectra 




Figure 6. Results from the 45 m and SMA observations of Galactic molecular gas along the sight line of MM J1545. (a) The 45-m 
integrated intensity map of (1-0) overlaid with the AzTEC 1.1-mm image (contours). The contours start at Icr (5.6 mjy beam~^) 

with an interval of Icr. The crosses mark the positions where (1—0) data were obtained, (b) The (1-0) and (1-0) spectra 

obtained with the 45 m telescope. The spectra are scaled by 0.3x for clarity. All of the spectra are shown in main-beam temperature 

scale Tmb S’IkI have a velocity resolution of 0.5 km s“^. With the 45-m beam (16''), (1-0) is clearly seen at the position of MM J1545, 

but the spectra at (0", 0") is consistent with the four adjacent spectra, suggesting that the emission does not come from MM J1545. 

No emission is found, (c) (d) The channel maps of (2—1) and (2-1) obtained with the SMA. The images are not cleaned. 

Crosses represent the position of the SMA 890 pm counterpart. The velocities in terms of the local standard of rest (LSR) are indicated at 
the top of each panel. The contours are drawn at (—4, —2, 2, 4, ...)xc7, where cr 0.1 Jy beam“^ with a velocity resolution of 0.55 km s“^. 
No compact emission is detected in (2-1) and (2-1). 


explained by strong gravitational lensin^ models (e.R., 
Negrello et al.l[2?np llWashi et al.ll2nill: IWardlow et alJ 

2013HVieira et al.ll2013HBussmann et al.ll2013D . In addi- 
tion, detailed modeling of the mm/submm source num¬ 
ber counts suggests that the excess of the counts at high 
flux densities (S'!.!!!!!!! > 10 mJy for example) are domi- 
nated by strongly leii sed SMGs as well as nearby galaxies 
(jNegrelio et al.ll20inH . 

The NIR sources, J1545B and J1541B, may be lensing 
objects. While the lack of multiband optical photome¬ 
try in J1545B makes it difficult to determine the lensing 
properties, the clear 2MASS/DSS detections of J1541B 
allow us to estimate the magnification factor and thus 
the intrinsic nature of MM-J1541. Hence, we focus on 
the lensing property of the MM-J1541-J1541B system 
first, and then discuss the MM-J1545-J1545B system. 

4.1.1. MM-JI 54 I 

To characterize the properties of J1541B, we per¬ 
form SED fits to the optical-to-NIR ^otometric data 
of J1541B using the Hyperz cod^^ developed by 
iBolzonella et al.l (200011. In SEP fits , we used SED tem- 
nlates of iBruzual fc Charl^ (1200,41) . Erom the SPIRE 
Av map, the visual extinction toward MM-J1541 is ap¬ 
proximately 1.5, but t he stray light from the Moon in the 
SPIRE observations (iRvgl et al.l 1201311 can bias against 
low visual extinction. We assume a conservative value 
of Ay = 1. Note that this assumption does not dra¬ 
matically affect the result because the total extinction is 
considered by combining the Galactic reddening and the 
intrinsic extinction, which eventually compensates the 
uncertainty in the Galactic extinction. Figure [7] shows 
the result of SED fits. The photometric redshift (photo- 

http://webast.ast.obs-mip.fr/hyperz/ 


z) of J1541B is Zh = 0.26to;i3 (x^ = 0.10, the error bar 
is from the 68% confidence interval). The optical to NIR 
SED is well described by a 33 Myr elliptical model with 
a stellar mass of 1 x 10^^ Mq and an intrinsic extinction 
of Ay = 0.10. The R-band decrement can be accounted 
for by the 4000 A break at z ^ 0.3, which is well within 
the photo-z range. The best-fit stellar age is young (33 
Myr) but maturer stellar SEDs can reasonably match the 
actual SED as well. 

Then we make use of the Faber-Jackson relation (FJR, 
iFaber fc Jacksonlll976ll to constrain the velocity disper¬ 
sion of the lensing galaxy, which allows us to model a 
singular isothermal sphere (SIS) as a lensing dark halo. 
FJR has originally been proposed as an empirical rela¬ 
tion between R-band absolute magnitude and velocity 
dispersion of galaxies. Similar relations a re now con¬ 
firmed in many filter bands o ut to the NIR (|Pahre et al.l 
II998I: iLa Barbera et al.ll201(Tl ). In our case, the redden¬ 
ing can be significant because of the foreground Galactic 
molecular cloud, so we use the NIR versi on of FJR pre¬ 
sented in Eq. 5 of iLa Barbera et'af] (|20I0D . The Einstein 
radius of a SIS mass distribution is expressed as 


= I" 154 X 


( _ 

V 200 kms“i / 




1 Gpc y \ 1 Gpc 




where is the velocity dispersion, R® and R^® are 
the angular diameter distances from the observer to the 
background source and from the foreground lens to the 
background source, respectively. In the redshift range of 
0.1 < zl < 0.7, which covers the 68% redshift confidence 
interval of J1541B, we find the velocity dispersion rang¬ 
ing 100 < (Jy < 250 km s“^ from the Ks magnitude and 
F JR. W hen we consider the background SMG at z = 3.0 
(S 13.2|) . then we find 0 e 0'.'3-2", which is smaller than 
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the actual separation between J1541B and the 24 /rm 
peak. We estimate the magnification factor /ig at the 
24 fim position (3'' apart from the len s centroid) u sing 
a gravitational leasing model GLAFIC (|Ogur lEoia), in 
which we employ a circular-symmetric SIS with the ve¬ 
locity dispersion cr^ derived above. For zl = 0.26 and 
zs = 3.0, we find 9e = 0"7 and the magnification fac¬ 
tor at the position of MM-J1541, /ig = 1.2. If we go to 
zl = 0.55 that is the edge of the 68 % confidence interval, 
then we have 0 e = 1”2 and /Xg = 1.4. To get the Ein¬ 
stein radius closer to the 3" separation angle, it would be 
necessary to double the angular diameter distance to the 
lens and/or increase the luminosity of the lens galaxy by 
a factor of ~ 4, which is very difficult to achieve within 
the uncertainties in the measured quantities. We have 
another likelihood peak at zl = ixt = 0 . 20 , 

see Figure 13 ), but this is unlikely because the inferred 
stellar mass is too high (1 x 10^^ Mq). Even if this would 
be the case, the magnification still remains a moderate 
value of fig = 2.0. Consequently, the magnification fac¬ 
tor of MM-J1541 is not likely as high as ~ 10 but rather 
moderate {fig ~ 1.2), suggesting that MM-J1541 might 
be an intrinsically hyper-luminous star-forming galaxy 
with a demagnified 1 . 1 -mm flux density of ^20 mJy or 
the intrinsic FIR luminosity of log(LFiR/T 0 ) ^ 13.7- 
13.9. 

4.1.2. MM-JI 545 

On the other hand, amplification for MM-J1545 at 
z ~ 4-5 may be larger although measurement of its mag¬ 
nification using existing data is difficult. The Fourier 
analysis of the SMA visibility data of MM-J1545 clearly 
shows its extended morphology (FWHM ^ l"-2", see 
Figure 13, suggesting that the source could be magni¬ 
fied by a foreground galaxy seen as J1545B. Unfortu¬ 
nately, we only have the limited number of photometric 
data points for J1545B, which likely suffer from relatively 
large Galactic reddening of Ay ~ 2.4. We find no ap¬ 
parent spectral break at Aobs > 1-2 fim, suggesting the 
lens redshift of zl < 2. So we assume several l ens r ed- 
shifts and a source redshift (zs = 4.8, from § 13.21) to 
estimate the Einstein radius in the same manner as MM- 
J1541. Note that a rest-frame ATs-band extinction of 0.2 
mag is used to correct the Galactic reddening. We find 
> 1”5 at Zl > 0.7, which is inconsistent with our 
SMA image showing the smaller separation (0'.'9) and no 
counter-images split by a gravitational lens. This sug¬ 
gest the lens redshift of z < 0.6. Actually, if the leasing 
galaxy is at zl = 0.5, we have the absolute magnitude 
in rest-frame ATg-band of Mk^ ~ —25 and ^ 160-170 
km s“^ from the FJR. This yields the Einstein radius of 
~ 0'.'9, consistent with the observed situation. 

Unfortunately, we cannot exactly predict the mag¬ 
nification factor of MM-J1545; the SMA 890-/im im¬ 
age exhibits neither multiple counter-images nor a large 
arc/ring, and the spatial extent of the source is un¬ 
known. Furthermore, the accurate redshifts of the lens 
and source are not available. All of the facts prevent 
us from precisely constructing a lens model. However, 
given that no lensed source with a magnification of 
/ig > 10 showing a single i mage with a l"-beam has 
been reported thus far (e.g., iBussmann et al.l[2013ll . it 
should be reasonable to assume fig ~ 10 as an up¬ 


per limit for the magnification factor of MM-J1545. In 
this case, the intrinsic FIR luminosity is still very high 
(4m > 1 X 1013 L©), even after correcting for mag¬ 
nification. Such a starburst galaxy in the hyperlumi- 
nous regime at z 4-5 is still rare compared to exist¬ 
ing studies le g. iRiechers et al.ll2010l: iWalter et al.l[20U^ 
iGombes et al.ll2012HVieira et al.ll2013l : iWeifi et al.ll2013H 

and is a unique laboratory to investigate properties of 
star-formation in the early Universe. 


o 



Wavelength (/xm) 

Figure 7. The optical to near-infrared spectral energy distribu¬ 
tion (SED) of J1541B. We fix a Galactic extinction of E(B — V) = 
0.3, which approximately corresponds to Ay = 1. The photometry 
from WISE (3.2 and 4.6 pm, open symbols) is not used because 
the WISE bands can be affected by emission from small grains 
and polycyclic aromatic hydrocarbons of J1541B. The stellar emis¬ 
sion from the background MM-J1541 can also contribute to the 
WISE photometry. The solid curve represents the best-fit SED at 
2 = 0.26^q^ 3 (68% confidence interval). The dashed curve shows 
the SED at the secondary minimum {z = 1.73^Q'g2)- The inset 
panel shows the reduced a function of redshift 2 :. 


4.2. The 1.1-mm Number Counts at the Brightest End 

The detections of ultra-bright sources allow us to con¬ 
strain the brightest end of the 1.1 mm number counts, 
which is complementary to the deep number counts 
at 5i.imm — 1-20 mJy obtained from our own 1.1- 
mm s urveys of SMGs (jHatsukade et al.l[201lHScott et al.l 
I20121) . The area where the la sensitivities of the Lupus- 
I AzTEC map are below 7 mJy beam“^ (typically 
5 mJy beam“^) is 3.65 deg^. We eliminate the high col¬ 
umn region where A y > 1 mag on the 2MASS extinction 
map (lDobashill2011h . which leaves 2.88 deg^. We detect 
three > 5 ct point sources over the 2.88 deg^ area (Tsuk- 
agoshi et al., in preparation). We care fully cross-identif y 
known starless and protostellar cores (iRvgl et al.ll2013ll . 
which leaves only two extragalactic sources, MM-J1545 
and MM-J1541. The inferred cumulative number counts 
-V(5'i.imm > 25 mJy) for Lupus-I are 0.69^045 deg“^ 
(the erro r is taken from the Icr confidence interval com¬ 
puted bv lGehreUll986ll . 

Furthermore, two additional brig ht sources (5'i ijnm = 
37.3 ± 0.7 mJy and 4 3.3 ± 8.4 mJy, llkarashi et ajJ^OlTl : 
iTakekoshi et al.l[2ni3[) have been found to date ove r the 
AzTEC survey fields (1.60 deg^, iScott et 3(ll2012[) and 
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Figure 8. The 1.1-mm cumu lative number counts obtained in the AzTEC surveys (filled circles for this wor k, diamonds for blank field 
surveys from lHcott et all i 20121). crosses for the AzTEC Cluster Enviro nmental Survey fromIScott et all II2Q12I) ). For co mparison, we plot 
the SPT 1.4-mm counts i Vieira^^t ajJ J2010l) . SCUBA S50-fi m counts l^Cnud^^^taL I 120081) . LABOCA 870-//m counts (jJohanssoi^^L 

S AL MA 870-//m OKarim ^^aLll^lSl) and 1.3-mm counts HHatsukade et aLII2013l) . We also show model predictions from lShimizi^^^ 
and IBethermin e^aLTH^cT^L the latter of which accounts for the strong-lensing effects. The 1.4-mm, 1.3-mm, 870-/im, and 850- 
fim counts are scaled to an equivalent 1.1-mm flux density using scaling factors of Si,imm/Si, 4 mm = 1-89, »S'i.imm/5'i.3mm = 1-41, 
Si.imm/Ss 70 (im = 0.54, and 5'i.imm/5'850^m = 0.51, respectively. 


the SMC peripheral field (1.21 deg^, iTakekoshi et al.l 
nnn), respectively. Both of them are also indicative 
of strongly-lensed magnification. Taking into account 
all of these four brightest extragalactic sources with 
■S'l.imm > 25 mJy, we estimate the cumulative number 
counts at > 25, 32.5, 40 mJy are 0.53;(!;^g, 

and 0.35;(!;23 deg respectively, where the la errors 
are again taken from iGehrelsI (|1986ll . 

Figure [8] shows the 1.1-mm cumulative number counts 
of the brightest AzTEC sources, as well as the deeper 
1.1-mm number counts obtained from blank fields 
toward the Great Observatories Origins Deep Sur¬ 
vey (GOODS) North and South fields, the Lockman 
Hole, the Gosmic Evolution Survey (COSMOS), the 
Subaru/XM M-Newton Deep Field, the Akari Deep 
Field-South (|Scott et al.l 120121) and from the AzTEC 
Cluste r Environmental S urvey (ACES) (IWilson et al.l 
I2008bt IScott et all 120121) . These counts intersect at 
‘S'l.imm 20 mJy, but the slope at S'l.imm ^ 20 mJy 
is shallower than the deeper part of the source counts. 
We also plot the 1.4-mm counts obtained by the 
South Pole Telescop e (SPT) survey (|Vieira et al.fl20ldl : 
IMocanu et al.l 120131) scaled to equivalent 1.1-mm flux 
densities using a s caling factor of S'i,imm/<5'i.4mm = 1-89 
(|Scott et al.ll2012D . The SPT counts do not include any 
nearby galaxies discovered by the Infrared Astronomical 
Satellite (IRAS) and those with synchrotron-dominated 
SEDs. Alternatively, these are thought to be dominated 
by strongly-lensed SMGs. The amplitude and slope at 
the brightest end of the 1.1-mm counts are consistent 


with those of the SPT counts given the uncertainty in 
the scaling facto r . Fur thermore, a model prediction by 
IBethermin et al.l (I2012D . where flux magnification due to 
strong gravitational leasing is accounted for, reproduces 
the 1.1-mm counts at S'l.imm > 25 mJy, providing more 
supporting evidence that the brightest sources found at 
1.1 mm are mostly attributed to the strong-lensing effect. 

4.3. The AzTEC-SPIRE Colors 

It is increasingly becoming clear that mm to long- 
submm observations are more likely to select higher red- 
shift SMGeHmnETRa^^hort-submm as predicted ear¬ 
lier (|Blain &: Longaii1ll993f) . and recent studies of ultra- 
bright sources selected in the mm to long submm ac¬ 
tually reveal a lot of SMGs at redshift z > 4 out to 
(|Vieira et al.ll20l4 iWeifi et al.ll20l4 iBoone et al.l 
I2013D . It has been also suggested that selecting red ob¬ 
jects in SPIRE bands whose SEDs are peaked at the 500- 
/rm band (i.e., objects that follows S'soo/im > > 

>S' 250 /im; so-called 500-/rm peakers) is also a useful way to 
pick up high-z candidate s even at z ^ 6 (jRiechers et al.l 
120131 : IDo well et ani2014D . MM-J1545 is estimated to be 
at z ~ 4-5 and is formally consistent with a 500-^m 
peaker, implying that 1.1-mm selection and (sub)mm 
color investigation are quite useful for isolating 2 > 4 
SMGs. 

In the left side of Figure[51 we plot the S'i,imm/<S' 350 /im 
and S'l.imm/'S'soo/im flux density ratios of MM-J1545 and 
MM-J1541. For comparison, we also show the same 
color-color plots of 48 AzTEG-selected sources with S/N 
> 4.0 from the AzTEG/ASTE GOODS-South survey 
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Figure 9. (Left) The S'l.imm/'S'sso^im versus 5i.imm/<S'500/im color-color diagram for MM-J1545 and MM-J1541 (large filled circle with 
error bars), along with 1.1-mm selected SMGs found in the GOODS-South field (dots). We also show predicted colors of template SEDs 
redshifted from 2 = 0 to 6 (colored in blue to orange); a gray-body with ^ = 1.8 and a dust te mperature of T^ust = 35 K (filled 
squares), Arp 220 (cros ses), M 82 (asterisks), SMM 42135—0201 (open circles. FSwinbank et al.ll20ig) . and a mean SMG (open squares, 
IMichalowski et al.ll201(ll . The SMGs of GOODS-South with spectroscopic redshifts are colored in the same way according to their redshift. 
The typical error bar of the GOODS-South flux density ratios is shown at the bottom right corner, which is mostly dominated by confusion 
noises in the SPI RE image. For refere nce, we also plot the flux density ratios for the most distant SMG known, HFLS3 at 2 : = 6.34 
(blue filled circle. [Riechers et al.ll2Ql^ . (Right) The Herschel/SPIKE S'35o^m/‘5250Aim versus Ssoonm/Ssso^m color-color diagram for 
MM-J1545 and MM-J1541. We also plot 1.1-mm selected SMGs found in the GOODS-South field with 350 and 500 fim detections. The 
symbols and colors represent the same as the left panel. The region surrounded by the straight lines that satisfies Ss^o^m /S250^im > 1 
and Ssoonm/Sszoum > 1 shows the criteria for the 500-/^m peakers, i.e., S250fim < S^so^m < S'sooMm- We also plot the region where 
5'250Mm < S350fim < SsoOfim/l-S, in which 2 : ~ 6 SMGs may fall. 


(|Scott et al.ll2010t iDownes et al.l[20T^ . twelve of which 
have robust counterparts with spectroscopic re dshifts 
ranging from z = 0.037 to 4.76 (lYun et al.ll2ni2[ l. The 
Herschel data are retrieved from the HSA and the Icr 
noise levels are ~ 0.5-0.6 mJy beam“^ in all the SPIRE 
bands. The detection thresholds are set to 2crfaint = 7.6, 
9.2, 10.4 mJy beam“^ at 250, 350, and 500 /rm, respec¬ 
tively, where CTfaint is an underlyin g confusion limit af - 
ter removing bright SPIRE sources (jNguven et al.ll201fll l. 
Forty AzTEC sources are detected at 500 /xm while 8 
AzTEC sources do not have a significant counterpart in 
any of the SPIRE bands. The redshift tracks of a modi¬ 
fied black-body with Tdust = 35 K and typical starburst 
galaxies are also overlaid. MM-J1545 is situated in a 
region of the plot where z ~ 4-5 galaxies are actually 
observed or are expected from redshift tracks of SEE) 
models. Similarly, the 1.1-mm to 350-/j,m color of MM- 
J1541 is consistent with those found in z Ri 3 galaxies. 

In the right panel of Figure [9l we show the 
*^ 350 /im/‘S' 250 /im and S*5oo/im/»5*250/im plots for 1VI1VI-J1545 
and MM-J1541, as well as the 40 AzTEC-selected sources 
in GOODS-South which are detected at 500 /xm. De¬ 
spite large uncertainties in flux density ratios, at least 16 
of the GOODS-South sources (40%) are consistent with 
500-/xm peakers that are detected at least at both 350 
and 500 /xm, and up to 21 sources (53%) may be 500- 
/xm peakers if we include sources only detected at 500 
/xm. At least 3 sources meet the criterion, S' 250 /im < 
Sssofim < ■ S.r^nni/.m/1-3, which is u sed to select z > 6 can¬ 
didates bv iRiechers et al.l (I2013I1 . The FIR-to-mm color 
of MM-J1545 is overall consistent with those of the 500- 


/xm peakers but slightly bluer than the z ^ 6 criterion, 
while the color of MM-J1541 is in good agreement with 
those of z ^ 3 sources, supporting the redshift estimates 
discussed in §0 

As demonstrated by these color-color plots, 1.1-mm 
selected sources that are very red in the (sub)mm will 
offer a unique opportunity to investigate how frequently 
massive starbursts are triggered in the z > 4 universe, 
which places constraints on galaxy formation models for 
massive dusty starbursts. Future follow-up studies us¬ 
ing the Atacama Large Millimeter/Submillimeter Array 
(ALMA) are needed to investigate this further. 

4.4. Is MM-JI 545 a First Hydrostatic Core? 

How to securely distinguish between SMGs and 
prestellar cores is always an issue in identifying SMGs 
behind Galactic molecular clouds, because flux densi¬ 
ties and FIR-to-submm colors of the brightest SMGs 
looks similar to those found in the low-luminosity end of 
prestellar cores. While MM-J1541 is well isolated from 
the Lupus-I molecular cloud [Ay — 1) and thus likely an 
extragalactic source, MM-J1545 is closer to the molec¬ 
ular cloud {Ay — 2.4, see also Figured]) and is worthy 
of assessing the possibility that MM-J1545 is a Galactic 
source. We hereafte r assume a dist ance to the Lupus-I 
cloud of D = 150 pc (|Comer6nll2008ll . where a 1" size cor¬ 
responds to the physical scale of 7.3 x 10“"* pc or 150 AU. 

The concordan ce scenario of low-mass star formation 
(|Shu et al.l I1987L for a review) begins with a collapse 
of a sub-pc-sized gravitationally-bound molecular core. 
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a so-called prestellar coreF^. Prestellar cores are ini¬ 
tially optically-thin to the thermal dust emission and 
isothermally collaps es, and thus it is represented as a 
cold (~ 10 K, e.sf.. iMarsh et alJ 12014( 1 and moderately 
dense (n(H 2 ) ~ 10^ cm“^) isothermal system. Most of 
them are found in regions where high H 2 col umn densi¬ 
ties o f fV(H 2 ) > 10^^ cm“^ are observed le.g. lBvsfl et al.l 
l201.3f) . The prestellar phase is followe d by protostellar 
phases (a.k.a., Cla sses 0, I, II, and III. iShu et al.l 1198^ 

I Andre et al.lll993ll . in which mid-IR and/or NIR emis¬ 
sion powered by a central protostar is always visible. A 
hypothetical small adiabatic region that is extremely- 
dense (n(H 2 ) ^ I0^°-10^^ cm“^) and compact 1- 
10^ AU), a so-called first hydrostatic core (lLarsonlll969ll . 
may occur at the center of a prestellar core in its fi¬ 
nal stage just before protostar formation, although the 
hrst core is expected to be short-lived (^ 10^ yr; e.g., 
ISaigQ &: Tomisak^l2011l: iTomida et al.ll20I3ll and the ob¬ 

servational nature is st i ll controversial (e.g . , lEnoch et al.l 
[Mfil: IChen i~aII[20Tnl: IPineda et al.l[2mil . 

MM-J1545 is a very unusual source if this would be 
associated with the Lupus-I cloud. The stark constraint 
on 24-/rm flux density (< 0.3 mJy), the absence of a 
bright compact NIR source at the position of the SMA 
source or extended reflection nebulosity (as seen in edge- 
on Herbig-Halo objects like HH30), the low (effective) 
dust temperature of Tdust ~ 7.1 K (§ 13.1.11) . and the 
non-detection of outflows virtually rule out the pro¬ 
tostellar Class-0 and advanced phases, suggesting that 
MM-J1545 might be in the prestellar phase and thus the 
mass should be dominated by dense molecular gas. C^®0 
J = 1-0 and 2-1 emission lines toward this source (i.e., 
z = 0) were, however, not detected with the SMA and the 
NRO 45 m telescope (FigurelHl), which places a constraint 
on molecular gas mass of Mlte < 0.005 Mq (§ 13.1.11) . 
In contrast, the 1.1 mm flux density of 44 mJy yields 
M(H 2 ) ^ 0.1 Mq if this is a starless core [Tdust = 7.1 K, 
Kd = 0.1 (A/250/xm)“^ g“^ cm^ (jHildebrandl 1 1 98.3[ ) with 
j3 = 1.4, and the dust-to-gas mass ratio of 100 are as¬ 
sumed]. This discrepancy between gas and dust masses 
indicates that MM-J1545 is not a prestellar core. 

Furthermore, the SMA continuum observations at 
890 /im and 1.3 mm revealed a very com pact source, 
compared to usual prestellar cores (§ EH). As shown 
in Figure El its visibility distribution does not show any 
evidence for an extended structure like a dense gas en¬ 
velope, which is at least ~0.01 pc and typically ~ 0.1 pc 
in size, correspond ing to an angular size of ^ 100" (e.g. 
lOnishi et al.ll200^ . In Figure El we show Fourier trans¬ 
form of 890-/im brightness di stribution predicted for a 
critical Bonnor-Ebert sphere ()Ebertlll955l : lBonnoil[l95fifl 
as a realization of a starless core in the Lupus-I cloud. 
The gas temperature and central H 2 density is assumed 
to be 7.1 K and 5 x 10^ cm“^, respectively. The total 
gas mass of 0.04 Mq is chosen so that the total 890-/im 
flux density matches the observed one. The overall vis¬ 
ibility slope and amplitude at higher spatial frequencies 
are, however, largely deviated from the prediction of the 

A compact dense molecular cloud core which lacks a central 
(proto)star, whether it is gravitationally bound or unbound, is re¬ 
ferred to as a starless core. Thereafter we use the term, prestellar 
core, to distinguish the advanced protostellar objects, and for sim¬ 
plicity we neglect whether it is bound or not. 


Bonnor-Ebert model. 

A possible Galactic inte rpretation of MM-J1545 might 
be a first hydrostatic core (lLarsodll969HMasunaga et al.l 

ll998HMasunaga fc Inutsuka]l2000H in its very early phase, 

because its compact nature and the very cold (T 10 K) 
and low-luminosity (Tboi 10 “^Tq) SED are consis- 
tent with those exp ected in a mostly-naked first core 
(|Tomida et al.l[2M^ . In this case, C^®0 emission line 
is optically thick and the non-detection of C^®0 is ex¬ 
plained if a significant fraction of gas mass, which is in¬ 
ferred from the dust SED (^ 0.1 Mq), might be stored 
behind the C^®0 photosphere. In Figure El we show the 
Fourier components of a 890-/im brightness distribution 
predicted for a first hydrostatic core, whic h is computed 
by radiation hydrodynamic si mulations (ITomida et al.l 
l2010t ISaigQ fc Tomisak^ 1201111 . The first hydrostatic 
core is produced by gravitational collapse of a 0.3 Mq 
Bonnor-Ebert sphere, and the mass and inclination an¬ 
gle are chosen so that the overall visibility amplitudes 
match the observed data. Consequently, the unresolved 
cusp-like structure appeared at \/v? -\-v'^ > 100 kA and 
the extended envelope predicted from the first core model 
are consistent with the actual SMA visibility amplitudes. 

Note that the adjacent compact NIR source, J1545B, 
and the observed JVLA 6 cm flux density at the SMA 
position cannot be explained consistently with the other 
data even in the first hydrostatic core models. We also 
note that first hydrostatic cores should be extremely rare 
because of the short lifetime; only one first core out of 
^100-1000 dense starless cores is expected, if compar¬ 
ing the lifetime with the dynamical time scale of star¬ 
less cores (^0.1-1 Myr). We therefore conclude that the 
z ~ 4-5 lensed SMG is the most likely and naturally ex¬ 
plained by the SED and geometry/spatial extent of the 
multi-wavelength counterpart to MM-J1545. 

5. CONCLUSIONS 

We report serendipitous detections of two 1.1-mm 
bright extragalactic sources, MM J154506.4—344318 
(MM-J1545) and MM J154132.7-350320 (MM-J1541), 
in the AzTEC/ASTE survey of the Lupus-I star-forming 
region. These sources are likely located at cosmologi¬ 
cal distances at z ~ 4-5 (MM-J1545) and z ~ 3 (MM- 
J1541). The main results are as follows: 

1. MM-J1545 is the brightest {Si. 1mm — 43.9 ± 
5.6 mJy) of ~ 1400 SMGs identified through 
the whole AzTEC 1.1-mm galaxy surveys. SMA 
(890 /rm and 1.3 mm) interferometry confirms the 
exact position, and photometry from VLA, ATCA, 
NMA, and Herschel, in addition to AzTEC/ASTE 
and SMA, constrains the SED well, which is 
in good agreement with a single gray-body with 
Td ust/(I + z) = 7.1 ± 0.3 K and ^ = 1.4 ± 0.1 
(§ 13.1.11) . The SED fits to the photometry at 
> 24 /im indicate a redshift of z ~ 3.4-5.6 (a com¬ 
bination of 68% confidence intervals of SED fits 
using the Arp 22(1 , average SMG, and SMM J2135 
templates) (§ 13.21) . This is also supported by the 
(sub)mm color analysis, in which we show that 
MM-J1545 is situated in a region where z ~ 4-5 
dusty galaxies are expected on t he S'sooitm/'S'i.imm- 
<S'350/xm/-5'i,imm diagram (§|4E1)- A faint NIR ob¬ 
ject, J1545B, is identified 0'.'9 east of the SMA 890- 

















































14 


Tamura et al. 


/im peak, which is likely a foreground lensing object 
that amplifies MM-J1545. Although it is difficult 
to constrain the SED of J1545B, the Einstein ra¬ 
dius would be ~ 0'.'9 if the lens redshift is zl ~ 0.5. 
Even if this is the case and the magnification factor 
might be as high as 10, the demagnified EIR lumi¬ 
nosity is still extreme (Lfir 10 ^^ Lq), suggesting 
that MM-J1545 is intrinsically a hyper-luminous 
galaxy at z ~ 4-5 f? 14.11) . 

2. MM-J1541 is identified as the fourth brightest 1.1- 
mm source (5'i.inim = 27.1 ± 5.0 mJy) of the 
whole AzTEC surveys. This object is also seen 
in the SPIRE 250-500 /rm and MIPS 24 /rm bands 
IS 13.1.21) . The inferred photometric redshift ranges 
z = 2.1-4.1 (a combination of 68 % confidence in¬ 
tervals of SED fits using the Arp 220, M 82, av¬ 
erage SMG, and SMM J2135 templates) (§ 13.2p . 
Again, this is also su ppor ted by the AzTEC- 
SPIRE color analysis (S 14.31) . An optical/NlR im¬ 
age (J1541B) is clearly offset from the 24-/xm cen¬ 
troid by ~ 3", and the SED of J1541B is consis¬ 
tent with a bright galaxy at Zl = 0 . 26 tQ; 23 J sug¬ 
gesting that J1541B gravitationally magnifies the 
background SMG, MM-J1541. Gravitational lens 
modeling using the Faber-Jackson relation and a 
singular isothermal sphere shows the magnifica¬ 
tion factor of MM-J1541 is moderate (/ig ~ 1.2), 
suggesting that MM-J1541 might be an intrin¬ 
sically hyper-luminous star-forming galaxy with 
log(LFiR/A 0 ) 13.7-13.9 (§0]). 

3. The brightest-end (S'l.imm > 25 mJy) of the 1.1- 
mm cumulative number counts is constrained by 
MM-J1545 and MM-J1541, in addition to another 
two sources from the literature; N(> 25mJy) = 
0-70 ^o; 34 deg“^ A^(>32.5mJy) = O.bStgJg deg”^, 
and N{> 40mJy) = 0.35;° g deg-^. The slope 
at S'l.imm ^ 20 mJy is shallower than the deeper 
part of the source counts obtained from general 
deep fields such as GOODS and GOSMOS. The 
amplitude and slope at the brightest-end is con¬ 
sistent not only with that properly scaled from 
the 1.4-mm counts obtained by the South Pole 
Telescope survey, which are thought to be dom¬ 
inated by strongly-lensed SMGs, but also with a 
model prediction where flux magnification due to 
strong gravitational lensing is accounted for. This 
suggests that a substantial fraction of S'l.imm > 
25 mJy sources may be gravitationally amplified. 

4. The overall SED from the optical to the radio and 
the spatial structure of MM-J1545 are explained 
neither by a local prestellar core nor a protostellar 
object, although it is found toward a relatively high 
H 2 column region of the local molecular cloud. A 
possible explanation for a Galactic object might be 
a first hydrostatic core with no envelope structure. 
Even in this case, however, neither the 6 -cm contin¬ 
uum emission nor the NIR compact object, J1545B, 
0'.'9 away from MM-J1545 can be explained by any 
models for a first hydrostatic core. On the other 
hand, a z ~ 4-5 SMG strongly lensed by J1545B 


naturally accounts for all of the observed proper¬ 
ties. Hence we conclude that MM-J1545 is not a 
first hydrostatic core but a galaxy at a cosmological 
distance. 

Unexpectedly, an extremely-bright SMG at z > 3 and 
a Galactic low-mass dense starless core (e.g., an ex¬ 
posed first hydrostatic core) could be similar in the mid- 
infrared to millimeter spectral energy distributions and 
spatial structures at least at > 1". This indicates that it 
is necessary to distinguish the two possibilities by means 
of broad band photometry from the optical to centime¬ 
ter, when a compact object is identified toward Galac¬ 
tic star-forming regions. (Sub)millimeter spectroscopy 
and/or sub-arcsec imaging of the object will help to de¬ 
termine the redshift and the presence of gravitational 
magnification, which will be able to carried out using 
ALMA. 
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